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The effect of an insulating barrier located at a distance a from a NS quantum point contact 
is analyzed in this work. The Bogoliubov de Gennes equations are solved for NINS junctions (S: 
anysotropic superconductor, I: insulator and N: normal metal), where the NIN region is a quantum 
wire. For a ^ 0, bound states and resonances in the differential conductance are predicted. These 
resonances depend on the symmetry of the pair potential, the strength of the insulating barrier and 
a. Our results show that in a NINS quantum point contact the number of resonances vary with the 
symmetry of the order parameter. This is to be contrasted with the results for the NINS junction, 
in which only the position of the resonances changes with the symmetry. 

PACS numbers: 74.20.Rp,74.50.+r,74.45.+c,81.07.Lk 



I. INTRODUCTION 



In high critical temperature superconductivity the 
symmetry of the pair potential is one the most studied 
aspects [ll 0] ■ Tunneling spectra depend strongly on this 
symmetry and therefore tunneling spectroscopy is a very 
sensitive probe for its study. In a d-symmetry and (110) 
orientation, for instance, the differential conductance has 
a peak at zero voltage, called zero-bias conductance peak 
(ZBCP) vifhich has been predicted theoretically by differ- 
ents works B, B, la l a B, Hi and observed experimentally 
lillSlllllllinii ■ The existence of the ZBCP is due 
to the formation of Andreev bound states at the Fermi 
level (zero energy states) near to the interface [l5llla.ll7| . 
These states appear due to the interference between scat- 
tering quasiparticles at the interface and the sign change 
of the pair potential. Studies of quantum point contacts 
in NIS junctions show that the ZBCP is removed by the 
quasiparticle diffractions at the point contact llsL [l9| . 
an aspect that has been shown experimentally |l4l |. Re- 
cently two quantum point contacts have been studied for 
the crossed Andreev reflection in d-wave superconductors 

On the other hand, in NINS jH ^ and NISN junc- 
tions 123, 13 , resonances in the differential conductance 
appear. In anisotropic superconductors, the resonance 
energies depend as well on the symmetry of the pair po- 
tential. In NINS junctions and (i2:y-symmetry, e.g. , the 
positions of these resonances are out of phase with re- 
spect to those predicted for isotropic superconductors |25| 
and in NISN junctions the conductance presents two 
kinds of resonances due to anisotropy of the pair poten- 
tial HI. 

In this paper, we analyze the differential conductance 



when quasiparticles are injected into a superconductor 
from a single-mode quantum wire, with an insulating 
barrier located at a distance a of the NS interface (NINS 
quantum point contact). We show that there exist bound 
states which cause resonances in the differential conduc- 
tance and that the number of these resonances depends 
on the symmetry of the order parameter. This is shown 
through the solution of the Bogoliubov-de Gennes equa- 
tion in NINS junctions, where NIN region is modeled by 
a wire of width W . In particular s and d - symmetries 
are considered. 



II. THE BOGOLIUBOV-DE GENNES 

EQUATION AND ITS SOLUTIONS IN NINS 

POINT CONTACTS 

The elementary excitations or quasiparticles in a 
superconductor are described by the Bogoliubov de 
Gennes (BdG) equations , which can be generalized for 
anisotropic superconductors ^^ . For a steady state these 
equations are 



Heiri)u{ri) + / dr2A(ri,r2)w(r2) = Eu{ri) 



-H;{ri)v{ri) + / dr2A(ri,r2)u(r2) = Ev{ri) 



(1) 



where He = —H^V^ /2m+V (r) — /x is an electronic hamil- 
tonian and fj. the chemical potential. A(ri, r2) is the pair 
potential, u(ri) and w(ri) are the wave function for the 
electron- and hole-like components of a quasiparticle, 



V'(r) 



u(r) 
v{r) 



(2) 
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It is suppossed that the quasiparticle moves on the x- 
y plane, the interfaces are normal to the a;-axis and the 
NIN region has a width W in the y direction, see fig. 



n The insulating barrier is modeled by a delta function, 
V{x) — UoS{x + a). The solutions of the BdG equations 
in the N/, N// and in the superconducting regions, are 
respectively, 
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FIG. 1: (a) The point contact NINS junction, the insulating 
barrier is located in x = —a and the NIN region is a single 
mode quantum wire with width W. (h) Scattering processes, 
the solid line and dashed line represent the electron and the 
hole-like components of a quasiparticle respectively. 
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where 



kf = ^kf ± 2mE/t9, fci = y'fcl - 7r2/W^2^ 
fc^(s) = 'Jkl±2mVL±{s)/h?. 

.i(y) = ^/^sin ^(y+Y 



^±(s) = v^'-iA±(s)r, 



(6) 

(7) 
(8) 



Mq (s) = 



E 



4{^) 



^±{s) 



E 



(9) 

The quasiparticles with k]^ and A:_ wavenumber move 

in the pair potential A+ and A_ respectively and given 

by 



A±(s) = A(±fc|i + sj) = A±e 



iipi 



with 



A(k) - /"d(ri -r2)e''"("-^-"-^)A(ri -r2) 



(10) 



All the evanescent modes have been neglected. This 
approximation is justified due to the fact that for tt < 
Wkp < 27r the narrow wire is a single mode and the 
energy of evanescent modes is well above the Fermi en- 
ergy 111 m. One finds A, B, C, D, E, F, c and d 
using boundary condictions in a; = —a and a; = 0. The 
electron-electron and electron-hole reflection coefficients 
are respectively 



Re 



, Rh 



2F. 



(11) 



where 



g ={1 + Z^) [{1 + Fif ^ F2F3] 

+ Z^ [(1 -Fif - F2F3] e-2«(fc+-fe-)« 

+ Z[1-F^ + F2F3] [Z{e 

+i{e 



2ik-a I —2ikj^a\ 



-2ikj^a 2ik-a\ 



(12) 



h ={Ff - F2F3 - I)[z2g2^(fe++fc_)a _ (^ _ j^)2] 

- Z(Z -I- i) [2Fi{e^'''+'' - e^'^-") 

+ (1 + F^ - i^2F3)(e2'^+'^ + e'*'=-°)] , (13) 



E, 



I 2 2 

^_^ dq ^^^\l^ cos^[7rg/2]/.(g), 

712^72, _ 1 7_^^ (i-g2)2 



(14) 
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(15) 



(16) 



III. DIFERENTIAL CONDUCTANCE 

Using the BTK model |23 , the normahzed differential 
conductance, Gj^, is calculated from 



Gf 



Gs ^ [(1 + Fpf 
Gn 



AZ^\ (1 - i?e + Rh 



4Fn 



(17) 



where Gjv the conductance when A = and a = 0. F^ 
is defined by H14() with fi = 1. For d-symmetry A± ~ 
Aq cos{2{9 ^ a)) , a is the angle between the (100) axis of 
the superconductor and the normal to the interface, and 
9 = sm^\s/kF). (cf. fig. nil.) 




FIG. 2: Diferenttal conductance for s -symmetry, (a) Different 
values of Z with a = 63W ; (b) different values of a with 
Z = 5. In both cases kpW = 1.7. 

Figures 13 and show the differential conductance for 
s and d2.2_j^2 -symmetries. When a = (NIS point con- 
tact), our results agree with ^^. For a ^ 0, subgap res- 
onances appear in the differential conductance and their 
number increases with a. When Z decreases, it can be 
seen that the number and position of the resonances re- 
main unchanged, only the peaks become broaden and for 
(i2,2_y2-symmetry the broad is greater for a fixed Z value. 

Figure 01 exhibits the differential conductance for a = 
7r/4 ((i2:-j^-symmetry). ZBCP does no apppear because 
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FIG. 3: Same as in fig. \^for q = (d^2_y2 -symmetry). 



its Andreev reflection is zero. In the last case the wave 
functions in the channel are a superposition of two plane 
waves with wave numbers ky = ±7t/W, each wave ex- 
periences a pair potential phase and tt respectiv- 
ley and therefore the Andreev reflection coefhcient a{6) 
{Rh{d) = \a,{S)\ ) for each wave are opposite, the waves 
of the reflected holes interfere destructively and the An- 
dreev reflections vanish. In relation with the dx-y - 
symmetry the number of resonances decreases compared 
whit the s and d^2_,y2— symmetries. Additionally when 
Z decreases the number of resonances is constant and 
the peak broadens and its position is smoothly displaced 
toward the right. 

The subgap resonances in the differential conductance 
are a direct consequence of the quasibound states formed 
inside the energy gap. The energies and lifetime of these 
bound states are given by the poles of the current trans- 
mission amplitude. Setting g = in Eg 1 121 one finds these 
poles. A complex energy, E = En -\- iEj, is required in 
order to solve this equation, where En is the position of 
resonance and fi/2 \Ei\ is the lifetime of the quasibond 
states. 

The resonance positions Eji for s or d2.2_„2-symmetries 
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FIG. 4: Same as in fig. \^for a — 7r/4 (dx-y -symmetry). 



are given by 



En == Eo{m: - (j)),n^ 1, 2, . 



(18) 



and for dxy symmetry are determined from 

En = Eo{2mT-(l)'), n=l,2,... (19) 



In these equations 



En 



E, 



I -If' 



(20) 



and <f> , (j)' are phases that depend on Z ^ a and E. There- 
fore the number of resonances with E < max{A) for s 
or d3.2„j,2 -symmetries are approximately twice the cor- 
responding number of a da^y-symmetry. This is due to 
the fact that in the case of a (i2;i/-symmetry the An- 
dreev reflection is zero. Thereby the quantization of the 
bound states occurs when the quasiparticlcs travel in a 
closed path a distance equal to 2a in the x direction and 
En DC 2n7r/2a. (The quasiparticle is transmitted as an 
electron in a; = — a and reflected as an electron in a; = 0.) 
In the case of s-or d^2_j,2- symmetries the quasiparticlcs 
complete a closed path when they travel a distance 4a 




FIG. 5; Diferential conductance for different values of a with 
Z = 5, a = 63W and kpW = 1.7. 



in X direction and En oc 2n7r/4a. (The quasiparticle is 
transmitted in a; = as an electron, reflected as a hole 
in a; = 0, reflected as a hole in x — —a and finally 
reflected as an electron in x = 0.) Therefore one has 
that in this case the number of the quasibound states is 
approximately twice the corresponding number of dx~y 
symmetry. 

In order to determine the lifetime of the quasibound 
states, a semiclassical argumentation will be used. The 
lifetime r is defined as the time that a quasiparticle in 
the INS region requires to "scape" toward the Nj or S 
regions. For the s or d,j.2^y2- symmetries the time that a 
quasiparticle needs for around trip is 



T ■ 



4a 



2M 



HkoFi/iTi 



If N is the number of closed trips, r is given by 



TiV, 



N is obtained from 



N 



[Re-hRh-eRh-hRe-e] — 1/e 



(21) 



(22) 



(23) 



where Re-h, Rh-e are the electron- hole and hole- 
electron refiection coefficients respectively for Z = 
(point contact NS) and RIe~e, RIh~h the electron- 
electron and hole-hole reflection coefficients respectively 
for an insulating barrier (IN). From equations (|21|l . H22() 
and H23|l the lifetime is obtained as 



M 



1 



Er. 



1 



— 7ln(Z2i?,_^/(l + Z2))^ 

IF 



(24) 



where we have used the fact that Re-h = Rh-e and 
RIe-e = Rlh-h = Z'^l (1 + 2'2). Equation ^ is sim- 
ilar to that found for NINS junction with s-symmetry 



[2l| . Similarly for the (ia;j,-symnietry the lifetime is given 

by 



hd 



EpJl-lp 



2ln(Z2i?,_,/(l + Z2))' 



(25) 



with Re-e the electron-electron reflection coefficient for 
Z = 0. For the case of s— symmetry , and E < |A| , 
Re~h = 1- Therefore the lifetime increases with Z and 
tends to infinity for Z » 1, while the resonance width, 
2 \Ei\ « h/r — > 0, as is observed in fig. El For da;2_,y2 — 
symmetry the quasiparticles transmission is finite for 
E < Aq due to the anisotropy of the pair potential, 
Re-h < 1, and the lifetime increases with Z but is fi- 
nite for Z >> 1. This is observed in the width of the 
resonances in fig. |3| For the d^y— symmetry the behav- 
ior of the lifetime and the width of the resonances are 
similar to the case of d^2_y2- symmetry, see Fig. 0] For 
all cases, with E > Aq, the refiection coefficients are al- 
ways less that one, the lifetimes decrease and the widths 
of the resonances increase. 

Figure shows Gr for different values of a, when a 
change from to 7r/4 , some peaks begin to decrease 
and vanish for a ~ 0.207r. This happens because the 
Andreev reflections decrease and the electron-electron re- 
flection increases. For a = 7r/4 the Andreev reflections 
are zero and one has the conductance for d^y- symme- 
try. Similarly the values of the energy of the resonances 
move toward the left as a increases due to a change of 
the phase 4> in the solution of the equation g — 0. 



IV. CONCLUSIONS 

Our results show that in NINS point contacts the 
differential conductance have resonances due to bound 
states. The number of resonances depends on the sym- 
metry of the order parameter, in contrast to a NINS 
junction. In the latter case only the position of the reso- 
nances changes with the symmetry. The number of res- 
onances with E < max{A) (subgap resonances) for s 
or d2,2_y2-symmetries is approximately twice the corre- 
sponding number of the dxy-symiaetry. When a change 
from to 7r/2 some peaks dissapear because to Andreev 
reflection vanishes 

In the case of s— symmetry, the lifetime of quasibound 
states increases with the insulating barrier strength and 
is infinite for Z » 1. In contrast, for a d-symmetry the 
lifetime increases with Z but is finite for Z >> 1. This 
occurs because the quasiparticles transmission is different 
of zero for E < Aq in contrast to the case of s-symmetry, 
where the transnnssion is zero for E < Aq. Therefore 
the lifetime of the resonances decreases in d-symmetries 
and the width of the resonances increases. The results 
obtained in this work can be used to find the symmetry 
of high temperature superconductors in experiments of 
the type carried out in references 'l(f\ and jlj. 
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